Electron spin coherence in Si/SiGe quantum wells 
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The mechanisms limiting the spin coherence time of electrons are of great importance for spin- 
tronics. We present electron spin resonance (ESR) and transport measurements of six different two 
dimensional electron gases in silicon/silicon- germanium (Si/SiGe 2DEGs). The spin decoherence 
time T 2 * is presented in conjunction with the 2DEG density n e and momentum scattering time r p 
as measured from transport experiments. A pronounced dependence of T 2 * on the orientation of 
the applied magnetic field with respect to 2DEG layer is found which is not consistent with that 
expected from any mechanism described in the literature. 



Spin decoherence mechanisms are of fundamental im- 
portance to spintronics. Silicon is an excellent model 
system for studies of decoherence, and electron spins in 
silicon have long coherence times, making them partic- 
ularly attractive for applications. The electron spin co- 
herence time T2 for phosphorus-bound donor electrons 
in isotopically pure 28 Si has been measured to be as long 
as 14 ms at 7K, and extrapolates to on the order of 60 
ms for an isolated spin 1]. Loss and DiVincenzo Q pro- 
posed the use of spins of single electrons as quantum bits, 
and Kane 3] has discussed the advantages of working 
in silicon. Vrijen and Yablonovitch et al. 4] have ex- 
tended that approach to include donor bound electrons 
in silicon-germanium heterostructures, and schemes have 
been proposed for electron spin-based quantum compu- 
tation in silicon-germanium electron quantum dots P, li| • 

There have been only a few studies examining the spin 
coherence time of silicon/silicon-germanium two dimen- 
sional electron gases (Si/SiGe 2DEGs) 0, H E3, OH- In 
this paper we provide a comprehensive treatment of six 
different samples, using transport measurements to ex- 
tract the electron density and scattering time and ESR 
to measure T 2 * and to provide an indication of the spin 
decoherence mechanism. A detailed study indicates that 
the dominant decoherence mechanism is strongly depen- 
dent on the orientation of the magnetic field, but it is 
inconsistent with the published mechanisms. 

The Si/SiGe heterostructures are grown by ultrahigh 
vacuum chemical vapor deposition at the University of 
Wisconsin - Madison and at IBM- Watson [n|. The 
2DEG sits near the top of a strained Si layer grown on a 
strain-relaxed Sii_. T Ge^ buffer layer, as shown in Figure 
1(a) of reference jjj. Above the 2DEG is a Sii-^Ge^ off- 
set layer, followed by a phosphorus-doped dopant layer, 
and then a Sii-^Ge^ spacer layer capped with Si at the 
surface. Table Q] contains details for each sample. 




FIG. 1: Electrons in the quantum well move in the presence 
of an electric field between the dopants and the well, and 
therefore feel an effective in-plane magnetic field. 



Hall measurements are performed on all samples. Hall 
bars are etch defined and Ohmic contacts are made to 
the 2DEG by Au/Sb metal evaporation and annealing at 
400 °C for 10 minutes. These data are used to extract the 
electron density and mobility, and from the mobility we 
derive the momentum relaxation time r p , an important 
parameter in spin relaxation via spin-orbit and related 
interactions. The parameters reported in Table [I] have 
been corrected for a small parallel conduction path using 
the method of Kane et al. |23| , a nd in each case this 
correction was smaller than 1% jl3| . 

Electron spin resonance data were acquired with a 
Bruker ESP300E X-band spectrometer, using an Oxford 
Instruments ESR900 continuous flow cryostat to main- 
tain a sample temperature of 4.2K. Magnetic field cali- 
bration and tracking was done with an ER035M NMR 
Gaussmeter. 

The ESR spectra for all samples were measured as a 
function of the orientation of the applied magnetic field, 
determined by the angle between the magnetic field 
and the growth direction of the sample (see Fig. [I}. Fig- 
ures Efa) and (c) are two-dimensional maps of the ESR 
intensity as a function of magnetic field and orientation 
angle for two selected samples. The peak-to-peak ESR 
linewidths AH PP were extracted by fitting the lineshapes 
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TABLE I: Sample parameters and measurements. The first section of the table contains growth parameters. The next section 
contains results from transport (Hall) measurements: in this section r p is calculated from the mobility fi by the expression 
r p = mtfi/e. The last three columns contain ESR measurements: T 2 * is derived from the peak width using Eq. (J (g = 2.00 for 
all samples), A(15°) is the anisotropy at = 15° for each sample, and b is the measured quadratic coefficient (from Eq. ([!}). 



to the derivative of a Lorentzian (see insets in Fig. EJb) 
and (d)). These linewidths show a pronounced depen- 
dence on the orientation angle (see Fig. EJb) and (d)). 
The minimal ESR linewidth (at = 0) and the observed 
linewidth anisotropy for all samples are summarized in 
Table HllH 

The ESR linewidth AH PP is directly related to the co- 
herence time T| through the expression [l^ 
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where g is the Lande g-factor and \±b is the Bohr mag- 
neton. It has been proposed that the orientation de- 
pendence of T 2 * (and thus of AH PP ) in similar 2DEG 
structures results from a D'yakonov-Perel spin relaxation 
mechanism due to fluctuating Rashba fields Elec- 
trons in the quantum well leave behind positive charge 
on their donors, setting up an electric field that is in ad- 
dition to any interface electric field. The electrons in the 
well on the Fermi surface move with the Fermi velocity in 
this electric field and therefore feel an effective magnetic 
field in the plane of the 2DEG (see Fig.QJ. This field is 
called the Rashba field Hrb- Scattering of the 2D elec- 
trons results in a fluctuating field Hrb which is always 
in the 2DEG plane (see Fig.GJ. Therefore, when the ex- 
ternal magnetic field Bo is applied perpendicular to the 
2DEG (0 = 0), the fluctuating Hbr are perpendicular 
to Bo. However, when Bo is tilted with respect to the 
2DEG (0 7^ 0), a component of the fluctuating field ap- 
pears along Bo driving decoherence more rapidly, which 
results in an orientational dependence of T 2 * and the ESR 
linewidth AH PP through Eq. QJ. In general, there may 
be other contributions to the linewidth as well, and the 
spin coherence time T 2 * can be written as 



T* 



T BR 



rp/ ' 
1 2 



where 1/T BR is the Rashba term contribution, and l/X^ 
includes all other contributions to the linewidth (such as 
inhomogeneous broadening or other decoherence mecha- 
nisms). 



Two groups have derived expressions for T BR in the 
limit uo c r p cos ^> 1. Both can be written as 
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where 77 = 1/2 for [g[ and 77 = 2 for [l^, a is the Rashba 
coefficient (defined by the Hamiltonian 7Y = a(crxk j p)-n, 
a are the Pauli spin matrices), k (fcp) is the electron 
(Fermi) wavevector, r p is the momentum relaxation time, 
is the angle the magnetic field makes relative to the 
sample growth direction, uo c = eB/m* is the cyclotron 
frequency, and uol = g^BH/ft is the spin precession (Lar- 
mor) frequency [3 . The limit uo { 
Eq. (|2J is valid only for small angles 0. 

If 1 / T RB is the dominant term in 1 / T 2 * , then Equation 
[21 can be normalized as follows to obtain an expression 
that does not depend on the value of the Rashba param- 
eter a, |2a| but only on the momentum scattering time 
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The orientational dependence of the normalized ESR 
linewidths for all samples and Eq. ((SJ (using each sam- 
ple's transport-measured t p ) are plotted in Fig. |31 The 
observed anisotropics at small angles differ substantially 
from those that Eq. JSJ) predicts. To see this quantita- 
tively, and since this expression only applies for small 0, 
we Taylor expand Eq. (|3J) to 



A(0) = 1 + b0 2 



(4) 



where the quadratic coefficient b is a measure of how 
quickly the anisotropy increases with angle 0. For each 
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FIG. 2: Orientation map of the ESR signal from (a) sample uw-031203 and (c) sample ibm-01 where the x-axis is the magnetic 
field, the shading scale is the peak intensity, and the y-axis is the orientation angle. Lorentzian-fit peak width (see inset) of (b) 
sample uw-031203 and (d) sample ibm-01 as a function of orientation angle. 



sample the quadratic coefficient b can be extracted from 
a parabolic fit (of Eq. (@J)) to the data. A plot of the 
quadratic coefficient b as a function of the parameter r p 
with the measured samples is given in Fig. 0] for both 
values of 77. For all six samples the quadratic coefficients 
b differ substantially from both theoretical predictions. 
Even more striking, the maximum quadratic coefficient b 
that Eq. Q can give (for any value of r p ) is 1.07 rad -2 , 
which is nearly an order of magnitude smaller than that 
observed for five of the six samples. 

As Fig. 0] shows for these samples, the semi-classically 
derived l/T 2 BjR does not account for the observed 1/T 2 *. 
A fully quantum mechanical derivation of the Rashba 
decoherence 1/T^ R may be necessary, or there may be 
other mechanisms contributing to the linewidth through 
the additional component l/T^. If the latter is the case, 
then 1/T 2 is necessarily orientationally-dependent, since 
an orientationally-independent term cannot change the 
functional form of the anisotropy, and so cannot affect 
b. That is, the observed discrepancy must be due to an 
orientationally-dependent effect. Since silicon possesses 
an inverison symmetry, orientationally dependent mech- 
anisms originating from the anti-symmetric term in the 
Hamiltonian introduced by Dresselhaus ^3 (see refer- 
ences 18] and 8]) should not contribute to the linewidth. 
It is conceivable that the data are showing some kind 
of orientationally dependent inhomogeneous broadening. 
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FIG. 3: Anisotropy (normalized peak width) as a function of 
orientation angle for all samples, and Eq. (0 (dashed 77 = 2, 
solid 77 = 1/2) for each using r p as reported in Table [I] 



One way to test this would be to do pulsed EPR experi- 
ments, measuring T 2 instead of T 2 * and removing the sen- 
sitivity to inhomogeneous broadening (such as the static 
dipole-dipole interactions with residual 29 Si). The inter- 
actions to residual 29 Si can also be eliminated through 
the use of isotopically pure 28 Si in the quantum well. 

In summary, we have performed a combination of ESR 
and transport measurements on six Si/SiGe 2DEGs and 
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FIG. 4: The quadratic coefficient b of Eq. (0 (dashed 77 = 2, 
solid 77 = 1/2) near the origin as a function of the parameter 
r p , including the measured sample 6's at transport r p 's (see 
Table |3. 



characterized the orientation dependence of the ESR 
linewidths. We observe an orient at ionally dependent spin 
decoherence with an anisotropy larger than any current 
theory predicts. 
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[25] In many of the ESR data sets (e.g., the inset of Fig.^b)) 
there is a small peak near 3341 Gauss (in the region of 
Lande g- factor g w 2.0) with no orient at ional dependence 
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[26] The presumed origin of the Rashba field in these samples 
is charge left behind in the dopant layer and the break- 
down of the effective mass approximation at the sharp 
quantum well interface. In general, any charge asymme- 
try can give rise to an electric field that would drive 
decoherence in the same manner. There are four differ- 
ent types of asymmetries: (a) bulk inversion asymm etry 
(BIA) due to the unit cell of the growth materials |20|, 
(b) structural inversion asymmetry (SI A) due to growth 
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sources of electric fields. In our case, with an asymmet- 
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from the dopant layer as in (b) and (local) fluctuations 
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